Non-thermal leptogenesis with strongly hierarchical right handed neutrinos 
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Assuming the Dirac-type neutrino masses mo are related to quark or charged lepton masses, 
neutrino oscillation data indicate that right handed neutrino masses are in general strongly hierar- 
chical. In particular, if mo is similar to the up-type quark masses, the mass of the lightest right 
handed neutrino M\ < 10 6 GeV. We show that non-thermal leptogenesis by inflaton decay can yield 
sufficient baryon asymmetry despite this constraint, and discuss how the asymmetry is correlated 
with the low energy neutrino masses and CP-violating phases. 
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I. INTRODUCTION 

An attractive mechanism for generating the observed 
baryon asymmetry of the universe (BAU) is baryogenesis 
via leptogenesis [1]. In the seesaw model [2], the out- 
of-equilibrium decays of right handed (RH) neutrinos to 
lepton and Higgs fields create lepton asymmetry, which is 
partially converted to baryon asymmetry by electroweak 
sphaleron processes [3] . 

The RH neutrinos can be generated thermally after 
inflation, if their masses are comparable to or below the 
reheat temperature T r . The thermal leptogenesis sce- 
nario has the nice feature that the final asymmetry is 
independent of initial conditions and inflaton couplings. 
However, it requires T r > 10 9 GeV to generate the BAU 
[4, 5], which is problematic in supcrsymmctric (SUSY) 
models due to the gravitino constraint [6]. Non-thermal 
leptogenesis by inflaton decay is an alternative scenario 
that can work with lower values of T r (> 10 6 GeV) [7- 
9], These bounds can be saturated with Mi ~ T r and 
Mi > T r for the thermal and non-thermal scenarios re- 
spectively, where Mi is the lightest RH neutrino mass. 

The seesaw relation 



to = mjjM 
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where is the Dirac-type neutrino mass matrix, re- 
lates the RH neutrino mass matrix M to the low energy 
neutrino mass matrix to, given in the basis where the 
charged lepton mass matrix and gauge interactions are 
diagonal by 



m = U, 
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Here d v = diag{mi,m2,m^) 
tonic mixing matrix 



and £7pMNS [10] is the lep- 
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Cij = cos 9ij , Sij = sin 9ij , 5 is the CP-violating Dirac 
phase and Kq = diag(e lai ^ 2 ,e la2 ^ 2 ,l) contains the two 
CP-violating Majorana phases. 

In Refs. 11-15, thermal leptogenesis was analyzed with 
the assumption that mo is related to the mass matri- 
ces of quarks or charged leptons, as typically realized in 
grand unified theories. In this case the Dirac masses are 
hierarchical, and the Dirac left-handed rotation in the 
basis where the charged lepton mass matrix is diagonal 
(the leptonic analogue of E/ckm) is expected to be nearly 
diagonal or similar to Uckm- We will hereafter refer to 
these two assumptions as quark-lepton symmetry. 

Hierarchical Dirac masses indicate strongly hierarchi- 
cal RH neutrino masses [13, 16], and the resulting BAU 
is suppressed due to the low value of M\. In partic- 
ular, Mi < 10 6 GeV if rriD is similar to the up-type 
quark masses. In this letter we point out that sufficient 
asymmetry can nevertheless be generated through non- 
thermal leptogenesis by inflaton decay. The inflaton is 
assumed to decay predominantly to the next-to-lightest 
RH neutrino. The asymmetry resulting from decays of 
this neutrino is partially washed out since Mi < T r . The 
final asymmetry depends on the asymmetry per neutrino 
decay as well as how strong the washout is. 

The plan of the paper is as follows: We first review the 
structure of seesaw parameters and estimate the asymme- 
try and the washout assuming quark-lepton symmetry. 
Numerical examples are provided in separate sections for 
normal and inverted hierarchical (or quasi-degenerate) 
light neutrino masses. We discuss how the BAU is cor- 
related with the CP-violating phases and conclude with 
a summary of results and some brief remarks on thermal 
leptogenesis. 



II. SEESAW PARAMETERS AND 
LEPTOGENESIS 

In the basis where the RH neutrino mass matrix is 
diagonal, the Dirac mass matrix can be written as 

m D = Uld D U R , (4) 

do = diag{rn,r)\, ttidi, tud^). Eq. (1) then takes the form 
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in 



Uld D Wd D Ul , 



(5) 
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where 



w = u R d- R l ul 



(6) 



is the inverse mass matrix of the RH neutrinos in the 
basis where m D = U L d D , and d^ = diag(M 1 , M 2 , M 3 ). 
From Eq. (5) one obtains 
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As mentioned in the introduction, we are assuming 
ttidi <C m^)2 "C tod3, and the Dirac left-handed rotation 
Ul ~ t^CKM ~ 1- Elements of m w m generally have 
a much milder hierarchy compared to the Dirac masses. 
The matrix W then has a simple hierarchical structure, 
and is diagonalizcd by [13] 
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= arg Mi . 



d 23 = rhe^m^r 

d 13 = m ee m^ T 

di2 = m ee m w 

K = ^a 5 (e-*/ 2 ,e-^ 2 / 2 ,e-^ 3 / 2 ), 

Here the phases of RH neutrinos 4>i are in 
to keep Mi real. The mass eigenvalues are 



Mi 



(10) 

The large neutrino mixings can originate from the seesaw, 
despite both Ul and Ur being nearly diagonal [16]. 

To estimate the BAU, suppose the inflaton predomi- 
nantly decays into the i-th family RH neutrino TVj. The 
comoving number density Yn is given by 
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Br < 1 is the branching ratio of the inflaton <fi to Ni, the 
factor 2 assumes <j) — * 27Vj, is the inflaton mass, and 
we have used the instantaneous decay approximation. A 
more accurate calculation shows Yn to be « 25% larger 
[17]. The asymmetry resulting from the decays of Ni 
(assuming it decays promptly [8]) is then 



Y A < 



2T r \ej\ V 
m 



(12) 



where A = (1/3)5 — L, is the lepton asymmetry pro- 
duced per decay of N i: and rj is a washout factor. In 
the simplest scenario, Mi ^> T r and there is no washout 
(i] = 1). On the other hand, if Mi < T r , part of the 
asymmetry will be washed out due to Ai mediated in- 
verse decays and AL = 1 scatterings. For M 2 < T r , N 2 
mediated processes contribute to the washout as well. 
The A asymmetry is multiplied by a conversion factor 
(C w 12/37 for SM and C w 10/31 for MSSM) to obtain 
the BAU resulting from sphaleron processes at equilib- 
rium above the electroweak scale [18]. 

For hierarchical RH neutrino masses as in Eq. (10), 

3aMim atm 



< 



16nv 2 



(13) 



where the parameter a = 1 for non-SUSY and a = 
2/sin 2 /3 for SUSY (tan/3 = (H u )/(H d )) and v = 174 
GeV [9, 19, 20]. Eqs. (12, 13) imply that if the inflaton 
decays into N x , the WMAP best fit Y B0 = 8.7 x lO" 11 
[21] for the BAU requires 



Mi 



> 



1 
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1.3xl0 6 GeV, (14) 



since m0 fC l r (in effect [22]). 

If the Dirac masses are related to the up-type quark 
masses, Eq. (10) indicates that Mi is too light to gen- 
erate the BAU. We will therefore assume > 2M 2 so 
that <j> predominantly decays into A 2 instead of N\. Us- 
ing Eqs. (4, 9), (rajjmcjij ~ mDimrjj with coefficients 
involving elements of m. It follows that for hierarchical 
Dirac masses the dominant contribution to the asymme- 
try from the decays of A 2 involves N 3 in the loop [23]: 
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In the expression e 2jCt , the label i = 1,2,3 refers to the 
RH neutrino, and a — e, fj,, r to the lepton flavor that 
it decays into. For quark-lcpton symmetry (that is, also 
assuming Ul ~ 1), the Dirac mass matrix has the form 

/ 0{m D1 ) < m D2 « m D3 \ 

m D ~ 0(m D1 ) 0(m D2 ) < m D3 _ ( 16 ) 

K <D{m D1 ) (D(m D2 ) 0(m D3 ) 

with coefficients involving elements of to, and the terms 
above the main diagonal proportional to non-diagonal 
elements of Ul- It follows that the dominant term in Eq. 
(15) is 
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(17) 
(18) 
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where ip < 1 is an effective phase that depends on 
d v , C^pmnSi dp and U^. (The phases fa in Ur can be 
calculated using Eqs. (5, 6) given the above masses, mix- 
ings and phases.) 

To estimate rji (the washout involving Ni), we define 
the washout parameters 



to* 



\m Dc 



Mi 



(19) 



to* w 1.08 x 1CT 3 cV for non-SUSY and to* w 
(sin 2 /3)1.58 x 10~ 3 eV for SUSY. Note that lepton flavors 
should be treated separately for an accurate calculation 
of the washout [14, 24-26]. Setting a = r, the washout is 
given in the instantaneous decay approximation by [5, 26] 



7] iiT w exp 



/ --* 3 /Cl(*0j(*)#i,T^T<l* 



(20) 



z = Mi/T, z = Mi/T r , JCi is a modified Bessel function 
of the second kind, and 



Yi a — —A a []Y\fj , 



(21) 



with I denoting the lepton doublet. The value of A TT de- 
pends on which interactions are in thermal equilibrium 
[24]. For MSSM, A TT = 19/30 between (l + tan 2 /3) x 10 5 
GeV and (1 + tan 2 /?) x 10 9 GeV [27]. For SM, A TT = 
344/537 and 390/589 below and above 10 9 GeV respec- 
tively [26]. 1 

The function j(z) takes AL — 1 scatterings into ac- 
count. We will not attempt a detailed calculation which 
involves finite temperature effects. Instead, we will 
use j(z) — 1 to define ?7 max which underestimates the 
washout, and 



3(z) = 



lC 2 {z) ( 9m 



Ki{z) \8ir 2 v 2 z 



+ 1 



to define ?7 m in which overestimates the washout [5]. 

It is also required that AL = 2 processes mediated 
by RH neutrinos are out of equilibrium. As discussed in 
Ref. 8, it is sufficient to have T r < (TO at m/"ii) 2 10 13 5 
GeV provided 



(a/8irv 2 )j: a \m Da2 \ 2 M 2 

V /27T 2 ^/45T2/TO P 



>1, 



(23) 



where mp w 2.4 x 10 18 GeV is the reduced Planck 
scale, and the relativistic degrees of freedom g* = 106.75 
(228.75) for SM (MSSM). Using ^>2,« ~ m Btm , this 
condition corresponds to M 2 > T r /5. 



1 A more accurate analysis around this temperature should take 
quantum oscillations into account [25]. 



For M 2 > T r , Tn 2 » and we can use the following 
simplified equations [4, 27, 28]: 



Z ^z~ = —z P *-lTz P «' 



(24) 



ZX ^ = *(Z-1)XY A + 2V *P* X ^ (25) 

dz z sHzm^ 

- \z*K, x {z)3{z)K^ T A Tr Y± T 

- i 7 ( 7 z) 3 /C 1 ( 7 z)j(7^)^l,r^rY AT . 



Here z = M 2 /T,j = M 1 /M 2 , and 



Z =l_ F ^ x — ( Pr + P<t> 



1/2 



(26) 



AHp. r \ p r 

with p r = (M 2 / z) 4 g*7r 2 /30. The equations are solved 
from Zi = M 2 /T max to z f > 7"\ and Y B w CYa t (z/)- 2 

III. RESULTS FOR NH SPECTRUM 

In this section we assume a normal hierarchical (NH) 
spectrum of light neutrino masses [m^ w TO a t m , to 2 w 
toq, toi <C to 2 ). To simplify the discussion we also set 
?7l = 1 and S13 — 0. In this limit the RH neutrino masses 
are given by [11, 13] 
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and with |c?i2 1 = \diz\ = s\ 2 m 2 m^/2 we obtain 



£2, 
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(22) Using Eqs. (12, 28), 
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(28) 

(29) 
(30) 



167TW 2 S 2 2 

For Mi <C T,. we can take z = in Eq. (20) to obtain 

(31) 



»)l,max = exp 
?7l,min = exp 



-^-K hT A, 



3?r 9m 2 
~8~87rV 



Yl^Tl 12 is the maximum temperature attained 
just after the inflaton starts oscillating, at time Hj 1 . We took 
Hi = for the numerical calculation, but the results are not 
sensitive to Hi as long as T m ax is at least a few times larger than 
T r (Hj » T r 2 /m P ). 
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FIG. 1: YB,max versus the reheat temperature T r , for do = 
d u , Ul = 1 and S13 = 0. The horizontal band corresponds 
to the WMAP range (8.7 ± 0.3) x lCT 11 [21], the solid and 
the dotted curves are calculated using Eqs. (24-26) and Eq. 
(20) respectively. Filled: SUSY, unfilled: non-SUSY; upper 
bounds: j(z) — 1, lower bounds: Eq. (22). 



r\2 can be estimated by using Eq. (20), and becomes 
significant for z$ < 10. It follows from Eqs. (29, 30) that 
the maximum asymmetry is obtained for T r sa M2/IO 
and ?Ti0 ss 20T r . Taking the reheating phase into account 
by solving the Boltzmann equations gives similar results. 

A numerical example is shown in Fig. 1 , where we have 
used Eq. (30) and set do = d u = diag{m Ul m c , m t ) with 
the values m u = 1.5 MeV, m c = 0.43 GeV, mt = 150 
GeV (taken from Ref. 29, for a renormalization scale 
of 10 9 GeV), for which M 1 sa 6 x 10 5 GeV and M 2 « 
6 x 10 9 GeV. 3 Assuming m\ -C m 2 , we can ignore the 
contributions to m^ a that involve m l5 and it follows from 
Eq. (29) that Y B oc m x . 

While the washout due to N 2 is severe when M 2 ^ T r , 
the washout due to Ni is rather mild, of order 0.1, for 
U L = 1 and S13 = 0. This follows from Eq. (20), with 
rfi2,T ~ "1-3/2 compared to rhi. T « c\ 2 m 2 /2 (in the limit 
nil -C ^2,3). For Ul ~ C^ckm, there are additional, 
order 9 c m 3 /m 2 contributions to m lr , where 9c is the 
lcptonic analogue of the Cabibbo angle. The result then 
depends on the CP- violating phases of C/pmns, but on 
average the washout gets stronger. 

We also take into account an effect due to off-diagonal 
elements of A a p [30] . Namely, in case of a strong washout 
related to a large rhi tT , part of the asymmetry can still 
survive if mi j(I or rh\^ e ~ m». 4 Typically mi, e is the 
smallest washout parameter. For an estimate we can ig- 
nore rh\^ which is of order m^ T1 and modify the Boltz- 



3 We take S12 = l/v3, S23 = l/v2 and sin/3 Ri 1 in the numerical 
calculations. We also take 7713 = 0.06 eV and m,2 = 0.011 eV, 
roughly approximating renormalization group effects by increas- 
ing the neutrino mass scale 20%. 

4 This is also true for the washout due to N2, and part of the 
asymmetry can survive for M2 J$ T r . However, the maximum 
BAU is still obtained for M 2 w l0T r . 
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FIG. 2: Histograms for YB, ma x with do = d u , mi = 0.2m2 
and Ul = Uckm, calculated for SUSY with j(z) = 1. Filled: 
S13 = 0, unfilled: S13 = 0.2. 



mann equations by adding 



1 



j{jz) 3 JC 1 (jz)j(jz)K hT A Te Y A 



(32) 



to Eq. (25), and including an analogous equation for Ya c 
(with r «-> e and £2 e ~ 0). The final asymmetry is then 
Y B ^C(Y A ^z f ) + Y Ae (z f )). 

To estimate the probability distribution of Ys imax for 
Ul = Uqkm, we numerically solved the Boltzmann 
equations 5000 times with uniformly distributed random 
phases of J/pmns- We define YB jmax by taking = 2M2 
as in Eq. (29), but here we take mi = 0.2m 2 to be specific 
and include (p. In addition, the asymmetry is maximized 
by varying T r for each run. Fig. 2 shows the results for 
S13 = and S13 = 0.2. The percentage of runs yielding 
YB,max > Ybo was 38% and 32% for S13 = and S13 = 0.2 
respectively. Including the A Te term significantly alters 
the low end of the probability distribution for YB )ma x> 
but the effect on these percentages is only a few points. 5 

For S13 = 0, the peak at Yb ~ 10~ 12 results from 
mi, e having a relatively small deviation (mi )e ~ [1 + 
O(0c) + 0(9 c m 3 /m 2 )}sl 2 m 2 ). For s 13 ^ 0, there are 
additional contributions to mi. e (as well as mi. r ), and the 
probability distribution becomes more dispersed. Using 
Eqs. (19, 4, 9, 10), J7ii >T ~ \m eT \ 2 /\m ee \, which from Eqs. 
(8, 2) is given by 



ci2Si2e la2 m 2 + I s 13 e 



-iS , Bc_ 
+ V2 



777.3 



(33) 



The terms including mi and evi are subdominant. As- 
suming 9q is similar to the Cabibbo angle, mi )T is mini- 
mized for 5 w 7r. yB !max also depends on £2, T , which from 



5 Note that in Eq. (4), there are generally Majorana phases on 
the left side of Ur as well. These phases enter rhi a , leading to 
O(0c) corrections for a = e,fi. (They can be rotated away in 
the limit Ul = 1.) This does not affect the results appreciably. 
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FIG. 3: Vs.max versus the Dirac phase S, with d_o = d u , 
mi = 0.2m2, t/t = (7ckm and S13 = 0.2, calculated for SUSY 
with j(z) — 1. 
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FIG. 4: Histograms for is, max with d D = d u , Ul = U ckm 
and s 13 = 0, calculated for SUSY with j(z) = 1. Filled: IH 
spectrum, unfilled: QD spectrum. 



Eq. (18) is maximized for 5 ~ 0. For random Majorana 
phases, Yb, max > Ybo is most likely at 6 » 7r due to 
the exponential dependence on mi |T , however it remains 
possible for all values of S (Fig. 3). 

The values of S13 and S will be probed by neutrino 
beam experiments within a decade for S13 > 0.05 [31]. 
The value of «2 can in principle be probed by neutrino- 
less double beta decay experiments. However, for normal 
hierarchy the effective Majorana mass Kto^)! = |m ee | ~ 
\s\ 2 e la2 m 2 + s 2 3 e~ 2lS m 3 \ is too small to detect using cur- 
rent techniques. 

IV. RESULTS FOR IH AND QD SPECTRA 

For inverted hierarchical (IH) spectrum of neutrino 
masses, m 3 <C rn\ < m 2 ~ m a t m and in the limit 
m 3 -> 0, t/ L -> 1, 



12' , 12 



2 to 2 



2 \c\ 2 e %ai m x 



sl 2 e la2 m 2 \ 



(34) 



Taking S12 = 1/V3, JWi.t ranges from 4m a t m /3 for 
«2 = Oil + 7T to mQ/36m^ tm for a 2 = ot\. (In- 
cluding the Cabibbo mixing, mi )T for a?2 = ai becomes 
w a tm^c/4; which is still < to*.) As a result, the asym- 
metry is suppressed by a factor ~ 10 8 for a 2 = a\ + ir, 
but Y B , max > Y B o is possible if ati « a 2 . 

The RH neutrino masses are given in this limit by 



Mi 



TO 



M 2 



2mf, 2 



M 3 



'D3 



TOatm 

for ai ~ a 2 . With |di2| 



Matm 2m 3 

|di3 1 = mim 2 /2 we obtain 



(35) 



IQlTV 2 



(36) 



similar to Eq. (28). 

The asymmetry can only survive if ati ~ a 2 for quasi- 
degenerate (QD) spectra of neutrino masses as well, since 



terms involving 7713 in toi jT are suppressed either by S13 
or 9c- The RH neutrino masses Mj ~ m 2 Di /fh where 
to is the QD neutrino mass scale. Assuming a% ~ a 2 , 
| £ 2,t| ~ (3a / 1 lQitv 2 )m 2 D2 is maximized for a 2 ~ 7T. On the 
other hand, TOi. r ~ Q 2 c fh for a.2 ~ it and it is minimized 
for ev 2 rs 0. The maximum asymmetry is determined by 
the interplay of these two factors. 

In the numerical examples we used the following neu- 
trino masses. IH spectrum: m\ = 0.059 eV, to 2 = 0.06 
eV, to 3 = to 2 /5. QD spectrum: toi = 0.1 eV, TO2 = 
0.1006 eV, to 3 = 0.117 eV. (Similar results are obtained 
for inverted hierarchical QD masses.) The resulting prob- 
ability distribution of Ye, max is displayed in Fig. 4. The 
percentage of runs yielding Yb, max > Ybo was 31% and 
18% for IH and QD spectra respectively, for Ul = Uckm 
and S13 = 0. Since rh\ , r cx to, Yb, max decreases as 
to is increased, with the percentage of runs yielding 
Y B , max > Y B0 decreasing to 7% (3%) for toi =0.2 (0.3) 
eV. These percentages increase a few points if Ul ~ 1, 
and decrease a few points if si 3 « 0.2. 

The effective Majorana mass, given by 



K TO /3/3>l 



TOi 



-TO 2 



(37) 



for both IH and QD spectra, is maximized by the condi- 
tion a 2 ~ ol\. As shown in Fig. 5, YB, m ax > Ybo requires 
I ("1/3/3) I > 0.04 eV for IH. 6 This range of \{mpp)\ can be 
probed within a decade [32] . 



V. CONCLUSION 

In this paper we considered non-thermal leptogenesis 
by inflaton decay under the assumption that the Dirac- 
type neutrino mass matrix tob is related to the up- type 



6 Note that since we took the neutrino mass scale 20% larger at 
the leptogenesis scale, we scaled |(mafl)| down 20% in the figure 
to correspond to low energy values. 
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\<m m >\ (eV) 

FIG. 5: YB,max versus the effective Majorana mass |(m^3)|, 
with do — d u and Ul = Uckm, calculated for SUSY with 
j(z) = 1. Black (red): IH spectrum with S13 = (si3 = 0.2). 
Blue (green): QD spectrum with S13 = (S13 = 0.2). 

quark (or charged lepton) mass matrix. Following the 
approach of Ref. 13, we did not make any specific as- 
sumptions on the textures of these matrices, but rather 
considered the general structure that follows from fitting 
to the low energy data. In this approach the RH neutrino 
masses are almost always strongly hierarchical (oc d 2 D ). 

This strong hierarchy and the rest of our analysis fol- 
lows from the light neutrino mass matrix having a less 
hierarchical structure compared to do, except for spe- 
cific values of S13 and the E/pmns phases which occur 
very rarely in a random scan. 7 On the other hand, rajj 
can be constrained further in particular SO(10) or other 
GUT models with flavor symmetries. It can then have 
structures different from Eq. (16), and fitting to the low 
energy data then yields those specific values as predic- 
tions. In such cases which are beyond the scope of this 
paper, the RH neutrinos can be less hierarchical and non- 
thermal leptogenesis with Mi ^> T r also becomes possi- 
ble. 

Assuming strongly hierarchical RH neutrinos with 
du ~ d u , the matter asymmetry created by the decays 
of N 2 is partially washed out since Mi < T r , but can 
still account for the BAU. For either NH or IH spectra 
of light neutrino masses, Yg = Ybo requires 

/^mWnnr«^ (j^V 10 -» , (38) 
V M 2 J V m 4> J \ m 2 J \m D2 J 

with each term on the left < 1 (min(mj) = mj for NH, 
min(mi) = 7773 for IH). The case do — d u corresponds to 
Mi ~ 6 x 10 9 GeV. 8 Eq. (38) can then be satisfied with 



Since nothing is currently known about the values of these 
phases, we assumed a flat probability distribution for numeri- 
cal analysis. Such a distribution also follows from 'anarchy' [33]. 
More precisely, M2 varies depending on the f/pMNS phases and 
is in the range 10 97 -10 10 - 2 GeV for NH with sia = 0. For NH 
with S13 = 0.2 or IH, it varies within an order (two orders) of 



TO0 = 10 10 -10 n GeV and T r ~ 10 s GeV. This value of 
T r can be consistent with the gravitino constraint, while 
TYi(f) in the above range is possible in small field or hybrid 
inflation models. On the other hand, for simplest large 
field inflation models > 2 x 10 13 GeV. The upcoming 
Planck satellite can discriminate these classes of models 

For NH spectrum of light neutrino masses and S13 > 
0.1, sufficient asymmetry is most likely to be obtained if 
the £/pmns Dirac phase 5 ~ it (assuming Ul « C/ckm)- 
For IH or QD spectra of light neutrino masses, sufficient 
asymmetry can only be obtained if the ?7pmns Majo- 
rana phases are approximately equal to each other, im- 
plying 1(771/3^)1 ~ ?7i a t m for IH spectrum and larger for 
QD spectrum. The asymmetry decreases as the QD neu- 
trino mass scale is increased, and if 1(771/3^)1 > 0.2 eV, the 
leptogenesis scenario discussed here is strongly disfavored 
assuming do ~ d u . 

If we relate the Dirac masses to masses of the charged 
lcptons di = diag{m e ,m^^m T ) instead, mo 2 ~ m-^tan/3 
and Eq. (27) yields M 2 ~ (tan 2 0)2 x 10 8 GeV. Pro- 
vided tan/3 is large, it is then easier to satisfy Eq. (38), 
especially for non-SUSY where there is no gravitino con- 
straint on T r . For large tan/3 it also becomes possible to 
generate the BAU with the inflaton decaying to Ni, as 
Mi ~ (tan 2 /3)5 x 10 4 GeV can satisfy Eq. (14). 

Thermal leptogenesis where the asymmetry is created 
by the decays of N 2 was discussed in Refs. 20, 30, 35 as 
well as Refs. 14, 15 which also relate mr> to the up-type 
quark masses. It is difficult to obtain sufficient asym- 
metry in this case. For 777,2, t ~ 77i a t m /2, the bounds are 
T r > 10 10 GeV and M 2 > 5 x 10 10 GeV assuming that e 2 
is given by Eq. (13) with Mi replaced by M 2 , and that 
there is negligible washout from Ni [5, 20]. However, 
for quark-lepton symmetry e 2 is suppressed by the light- 
est neutrino mass, and the phase values that maximize 
it do not coincide with those that suppress the washout. 
We therefore expect these bounds to be at least a few 
times larger. Similar conclusions are reached in Refs. 
15. Notwithstanding the high T r , the value of M 2 would 
then not be compatible with the assumption dp ~ d u , 
although it may be compatible with do ~ di tan (3 for 
large tan/3. 
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